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Defect production due to quenching through a multicritical point 
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We study the generation of defects when a quantum spin system is quenched through a multi- 
critical point by changing a parameter of the Hamiltonian as t/r, where r is the characteristic time 
scale of quenching. We argue that when a quantum system is quenched across a multicritical point, 
the density of defects (n) in the final state is not necessarily given by the Kibble-Zurek scaling form 
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where d is the spatial dimension, and v and z are respectively the correlation 



length and dynamical exponent associated with the quantum critical point. We propose a general- 
ized scaling form of the defect density given by n ~ i/t-'*/(2z2)^ where the exponent Z2 determines 
the behavior of the off-diagonal term of the 2x2 Landau-Zener matrix at the multicritical point. 
This scaling is valid not only at a multicritical point but also at an ordinary critical point. 
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INTRODUCTION 



The zero temperature quantum phase transitions oc- 
curring in quantum many body systems has been a chal- 
lenging area of research for the past few yearsA^. The dy- 
namics taking place in such systems by varying a param- 
eter in the Hamiltonian in a definite fashion have come 
to the forefront only recently^iii^i^ii. In this paper, we 
focus on the density of defects generated when a system, 
prepared in its ground state, is adiabatically quenched at 
a uniform ratei^>^^>2iioaiii2a3a4^a6^a8^^^. These 
works have their roots originating from the study of phase 
transitions in the early universe^^ which was extended to 
second order phase transitions^^ and later to quantum 
spin chains^. The diverging relaxation time associated 
with a quantum critical point results in the failure of 
the system to follow its instantaneous ground state; this 
eventually leads to the generation of defects in the final 
state. When a parameter of the quantum Hamiltonian is 
varied as i/r, where t is the characteristic time scale of 
the quenching, the Kibble-Zurek (KZ) argument^ pre- 
dicts a density of defects in the final state that scales as 
-^ I ^du / {zv+i) jj-^ ^j^g limit r — > oo. Here v and z denote the 
correlation length and dynamical exponents, respectively, 
characterizing the associated quantum phase transition 
of the d-dimensional quantum system. The KZ predic- 
tion has been verified for various exactly solvable spin 
models when quenched across a critical point^i^i^i^. Var- 
ious generalizations of the KZ scaling form have also been 
proposed for quenching through a gapless phase or along 
a gapless lin e^°d^'^° . Experimental verification of the 
dynamics of such systems can be realized by the trapped 
ultra cold atoms in optical lattices; for a review see Ref. 



The generation of defects during the adiabatic quench- 
ing dynamics of a one-dimensional spin-1/2 XY chain 
across a quantum critical point was studied in refla- The 



Hamiltonian of the system is2^ 



where the cr's are Pauli spin matrices satisfying the usual 
commutation relations. The strength of the transverse 
field is denoted by h, and J^ and Jy are the strength of 
the interactions in the x and y directions, respectively. 
The phase digram of the above model is shown in Fig. 1 . 
It was observed that when the transverse field h is 
varied as ft. = i/r, the system crosses the Ising transi- 
tion lines as shown in the figure, and the defect density 
scales as^ l/v^- This is in agreement with the KZ pre- 
diction since the values of the critical exponents associ- 
ated with the Ising transition are given hy z = u = 1. 
On the other hand, if the interaction in the a:-direction 
{Jx) is quenched in a similar fashion keeping h and Jy 
fixed^, the defect density is again found to scale as 1/a/t 
though the magnitude depends upon the values of Jy and 
h. If /i < 2Jy, the system crosses the anisotropic critical 
line {Jx — Jy) in addition to the Ising transition lines 
mentioned above, and hence the magnitude of the de- 
fects is increased. However, it was observed that if Jx is 
quenched keeping h — iJy, the system crosses the mul- 
ticritical point at Jx = Jy and h = 2Jy, where the Ising 
and anisotropic transition lines meet. The density of de- 
fects in the final state generated in a passage through 
the above multicritical point shows a slower decay with 
T given as l/r^/^. Since the critical exponents associ- 
ated with this multicritical point are given by i' = 1/2 
and z = 2, the above scaling relation does not follow from 
the KZ scaling relation \ I t'^'^ / 'y'^'^+'^) . It is this observa- 
tion which motivated us to look for a generalized scaling 
relation valid even for a multicritical point. It should 
be noted here that this is the first attempt to provide a 
generalized scaling relation for defect density when the 
system is quenched linearly through a multicritical point 
which has also been extended to the nonlinear case in a 
recent work^S . 
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FIG. 1: The phase diagram of the anisotropic XY model in 
a transverse field in the h/{Jx + Jy) — 7 plane, where 
7 = {Jx — Jy)/{Jx + Jy)- The vertical bold lines given by 
h/{Jx + Jy) = ±1 denote the Ising transitions. The system is 
also gapless on the horizontal bold line 7 = for \h\ < Jx + Jy 
FMx (FMy) is a long-range ordered phase with ferromagnetic 
ordering in the x (y) direction. The thick dashed line marks 
the boundary between the commensurate and incommensu- 
rate ferromagnetic phases. The thin dotted lines indicate the 
adiabatic and impulse regions when the field h is quenched 
from —00 to cxD. The two points with coordinates 7 = and 
h/ (Jx + Jy) = ±1 denoted by A and B are multicritical points. 



The paper is organized as follows. In Sec. II, we derive 
the general form for the scaling of defects and apply it in 
two models. Sec. Ill consists of concluding remarks. 



II. GENERAL SCALING 

To propose a general scaling scheme valid even 
for a multicritical point using the Landau-Zener non- 
adiabatic transition probability— i^, let us consider a d- 
dimensional model Hamiltonian of the form 

7f = X^ ^t(fc)[ (A(i) + fe(fc) ) a^ 

k 

+ A(fc) a+ + A*{k) a- ] V(fc), (2) 

where a^ — {a^ ± icr"^), b{k) and A(fc) are model depen- 
dent functions, and ip{k) denotes the fermionic operators 
(^i(fc), 'ip2{k)). The above Hamiltonian can represent, for 
example, a one-dimensional transverse Ising or XY spin 
chain2^, or an extended Kitaev model in d = 2 written 
in terms of Jordan- Wigner fermions^^. We assume that 
the parameter X{t) varies linearly as t/r and vanishes 
at the quantum critical point at t = 0, so that the sys- 
tem crosses a gapless point at t — for the wave vector 



k = ko. Without loss of generality, we set |fco| =0. The 
parameters b{k) and A(fc) are assumed to vanish at the 
quantum critical point in a power-law fashion given by 



b{k) - |fc|"i and A(fc) - \k\' 



(3) 



Many of the models described by Eq. ^ exhibit a quan- 
tum phase transition with the exponents associated with 
the quantum critical point being v = z = Z2 = 1- We 
shall however explore the more general case below. 

The Schrodinger equation describing the time evo- 
lution of the system when A is quenched is given by 
idip/dt — Hij} (where we set Planck's constant h ^ 1). 
Using the Hamiltonian in Eq. ^, we can write 



b{k) Mk) + A(fc) Mk), 



= - - + ^(^) Mf") + ^*(^) ^i(^)- (4) 



' dt 
' dt 



One can now remove 6(fc) from the above equations 
by redefining t/r + b{k) — > t\ thus the exponent Zi 
defined in Eq. ([3]) does not play any role in the 
following calculations. Defining a new set of vari- 
ables il^i{k) = ■ipi{k)exp{i J dt' t'/r) and ?A2(fc) = 
ip2{k) exp{~i J dt' t'/r), we arrive at a time evolution 
equation for 'ipi{k) given by 



dt^ 



2i- 



.t d 
"di 



Further rescaling t 



\A{krJMk) 
> tr^/^ leads to 



0. 



dt^ ~ dt 



+ \A{k)\\]Mk) = 0. 



(5) 



(6) 



If the system is prepared in its ground state at the be- 
ginning of the quenching, i.e., ipiik) = 1 at t = —00, the 
above equation suggests that the probability of the non- 
adiabatic transition, pr — limt^+00 IV^iC^)!^; must have 
a functional dependence on |A(fc)pT of the form 

p^, = /(|A(fc)|V). (7) 

The analytical form of the function / is given by the 
general Landau-Zener formula^^'^-^. The defect density 
in the final state is therefore given by 



d'^k 
(2^ 



/(|A(fc)|V) ^ 



d'^k 



H !{\kr'r). (8) 



The scaling k -^ k'^^^T finally leads to a scaling of the 
defect density given by 



1/r 



d/{2z2) 



(9) 



We shall recall the example of the quenching dynamics 
of the transverse XY spin chain when the field or the 
interaction is quenched^>^. When the system is quenched 
across the Ising or anisotropic critical line by linearly 



changing h or Jx as t/r, A(fc) vanishes at the critical 
point as A(fc) ~ \k\ yielding Z2 — z — 1; hence the gen- 
eralized scaling form given in Eq. ([9]) matches with the 
Kibble-Zurek prediction with v = z = 1. On the other 
hand, when the system is swept across the multicritical 
point {Jx = Jy,h = 2Jy) by quenching the interaction 
Jx = t/r with h — 2Jy, the equivalent 2x2 Hamiltonian 
matrix of the Jordan- Wigner fermions in an appropriate 
basis can be written as^ 

Jx + Jy (cos 2fc -I- 2 cos k) Jy (sin 2fc -I- 2 sin k) 

Jy{sm2k + 2sinfc) ~Jx — Jy {cos2k + 2cosfc) 

The corresponding Schrodinger equations are 



' dt 



' dt 



+ Jj^(cos2fc + 2cosA:) ) 'ipi{k) 
Jy{sm2k + 2sinfc) V'2(fc), 
Jy (sin 2fc + 2 sin k) ipi{k) 



h Ji,(c0s2fc -I- 2c0s/c) ) "02 (^) 



(10) 



At the quantum critical point Jx — Jy, the diagonal term 



6(fc) — Jy {cos2k + 2cosfc) goes as — Jy 



Jy\TT 



fc|2 



near fc = tt. Hence the dynamical exponent is given by 
z = zi = 2 at this multicritical point. Note that in 
this example, the critical point is not crossed at i = 0; 
however, one can shift the time so that b (k) ~ JTr— fc|^\ 
which would ensure that the quantum critical point is 
crossed at i = 0. 

On the other hand, the off-diagonal term A(A;) = 
Jy(sin2/c-|-2sinfc) — JTr— fcp leads to the density of defect 
scaling as 1/t^/^; this is in agreement with the general- 
ized scaling relation proposed in Eq. ^ with Z2 = 3. 
Fig. 2 shows the numerical integration of Eq. pU|) which 
confirms the defect scaling exponent of —1/6. 

Finally, let us comment on the dynamics of an exactly 
solvable transverse Ising model with an additional three- 
spin interaction which is also quenched through a multi- 
critical point^ by varying the transverse field as /i = t/r. 
It has been observed that the defect density n scales as 
1/t^/^ which again does not support the KZ scaling form. 
The three-spin interacting Hamiltonian is given by'^*' 



^— Je-^ 



[h + J3<-i<+i] 
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FIG. 2: n vs T obtained by numerically solving Eq. (|10p at 
the multicritical point with h = 10 and Jy = 5. The line has 
a slope of —0.16. 



phase diagram of this model is at /i = — 1 and J3 = 1/2. 
We observe that the off-diagonal term sinfc — J3 sin2fc 
scales as j/cp at the multicritical point; therefore the de- 
fect density scales as 1/t^' ^ as expected from the general 
scaling relation proposed here. The importance of the 
multicritical point has also been observed in nonlinear 
quenching of different models^S. 



III. CONCLUSIONS 

We have shown that the density of defects n produced 
when a system is quenched through a multicritical point 
does not follow the KZ scahng relation l/r''''/(^'^+^^. We 
have then proved a new scaling form which is not only 
valid at an ordinary quantum critical point but is also 
valid at a multicritical point. We argue that for a sys- 
tem which is swept across a multicritical point in a phase 
diagram, it is the exponent Z2 defined above which ap- 
pears in the scaling of the defect density given in Eq. ^ . 
However, for a passage through an ordinary critical point 
in many models, Z2 = z = 1, and Eq. ([9]) reproduces the 
conventional KZ scaling form with ly = z = I. 



We shall henceforth set J^ = 1- The equivalent 2x2 
Hamiltonian matrix of the Jordan- Wigner fermions in 
the momentum representation takes the form 



h(t) + cos k — J3 cos 2k 
— i(sinfc — J3sin2fc) 



z(sin k — J3 sin 2k) 
-{h{t) -t- cos k — J3 cos 2k) 



It may be noted that by virtue of a duality transfor- 
mation, this model can be mapped to a transverse XY 
model with competing interactions for the x and y com- 
ponents of the spin^. The multicritical point in the 



Acknowledgments 



AD acknowledges Subir Sachdev for a very interesting 
discussion. AD and UD thank G. E. Santoro for stimu- 
lating discussions and the hospitality of SISSA, Trieste, 
Italy where part of this work was carried out. AD and DS 
also thank Krishnendu Sengupta for his comments. DS 
thanks DST, India for financial support under Project 
No. SR/S2/CMP-27/2006. 



* Electronic address: 

^ Electronic address: 

^ Electronic address: 

^ Electronic address: 



|udiva@iitk.ac.inj 
|victor@iitk.ac.in] 
dutta@iitk.ac.in 



diptiman@cts . iis c.ernet.in| 
S. Sachdev, Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambridge, 1999). 

B. K. Chakrabarti, A. Dutta, and P. Sen, Quantum Ising 
Phases and Transitions in Transverse Ising Models, m41 
(Springer- Verlag, Berlin, 1996). 

T. Kadowaki and H. Nishimori, Phys. Rev. E 58, 5355 
(1998); K. Sengupta, S. Powell, and S. Sachdev, Phys. Rev. 
A 69, 053616 (2004); P. Calabrese and J. Cardy, J. Stat. 
Mech: Theory Expt P04010 (2005), and Phys. Rev. Lett. 
96, 136801 (2006); A. Das, K. Sengupta, D. Sen, and B. 
K. Chakrabarti, Phys. Rev. B 74, 144423 (2006). 
W. H. Zurek, U. Dorner, and P. ZoUer, Phys. Rev. Lett. 95, 
105701 (2005); J. Dziarmaga, Phys. Rev. Lett. 95, 245701 
(2005); B. Damski, Phys. Rev. Lett. 95, 035701 (2005); B. 
Damski and W. H. Zurek, Phys. Rev. A 73, 063405 (2006) 
A. Polkovnikov, Phys. Rev. B 72, 161201(R) (2005). 
R. W. Cherng and L. S. Levitov, Phys. Rev. A 73, 043614 
(2006). 

R. Schiitzhold, M. Uhlmann, Y. Xu, and U. R. Fischer, 
Phys. Rev. Lett. 97, 200601 (2006); C. KoUath, A. M. 
Liiuchh, and E. Altman, Phys. Rev. Lett. 98, 180601 
(2007); S. R. Manmana, S. Wessel, R. M. Noack, and A. 
Muramatsu, Phys. Rev. Lett. 98, 210405 (2007); M. Eck- 
stein and M. KoUar, Phys. Rev. Lett. 100, 120404 (2008). 
V. Mukherjee, U. Divakaran, A. Dutta, and D. Sen, Phys. 
Rev. B 76, 174303 (2007). 

U. Divakaran and A. Dutta, J. Stat. Mech: Theory and 
Experiment PllOOl (2007). 

D. Sen, K. Sengupta, and S. Mondal, Phys. Rev. Lett. 101, 
016806 (2008); R. Barankov and A. Polkovnikov, Phys. 
Rev. Lett. 101, 076801 (2008). 

J. Dziarmaga, Phys. Rev. B 74, 064416 (2006); T. Caneva, 
R. Fazio, and G. E. Santoro, Phys. Rev. B 76, 144427 
(2007). 

F. M. Cucchietti, B. Damski, J. Dziarmaga, and W. H. 
Zurek, Phys. Rev. A 75, 023603 (2007). 
J. Dziarmaga, J. Meisner, and W. H. Zurek, Phys. Rev. 



Lett. 101, 115701 (2008). 
^'^ C. De Grandi, R. A. Barankov, and A. Polkovnikov, Phys. 

Rev. Lett. 101, 230402 (2008). 
^^ A. Polkovnikov and V. Gritsev, Nature Physics 4, 477 

(2008). 
^^ K. Sengupta, D. Sen, and S. Mondal, Phys. Rev. Lett. 100, 

077204 (2008); S. Mondal, D. Sen, and K. Sengupta, Phys. 

Rev. B 78, 045101 (2008). 
^'' F. Pellegrini, S. Montangero, G. E. Santoro, and R. Fazio, 

Phys. Rev. B 77, 140404(R) (2008). 
^^ D. Patane, A. Silvia, L. Amico, R. Fazio, and G. Santoro, 

Phys. Rev. Lett. 101, 175701 (2008). 
^^ V. Mukherjee, A. Dutta, and D. Sen, Phys. Rev. B 77, 

214427 (2008). 
^° U. Divakaran, A. Dutta, and D. Sen, Phys. Rev. B 78 

144301 (2008). 
^^ T. Caneva, R. Fazio, and G. E. Santoro, Phys. Rev. B 78, 

104426 (2008), and Phys. Rev. B 76, 144427 (2007). 
22 T. W. B. Kibble, J. Phys. A 9, 1387 (1976), Phys. Rep. 

67, 183 (1980). 
2=* W. H. Zurek, Nature (London) 317, 505 (1985), and Phys. 

Rep. 276, 177 (1996). 
2'' I. Bloch, J. Dahbard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
25 E. Lieb, T. Schuhz, and D. Mattis, Ann. Phys. (NY) 16, 

407 (1961); E. Barouch and B. M. McCoy, Phys. Rev. A 3, 

786 (1971); J. E. Bunder and R. H. McKenzie, Phys. Rev. 

B 60, 344 (1999). 
2"^ C. Zener, Proc. Roy. Soc. London Ser A 137, 696 (1932); 

L. D. Landau and E. M. Lifshitz, Quantum Mechanics: 

Non-relattvistic Theory, 2nd ed. (Pergamon Press, Oxford, 

1965). 
2^^ S. Suzuki and M. Okada, in Quantum Annealing and Re- 
lated Optimization Methods, Ed. by A. Das and B. K. 

Chakrabarti (Springer- Verlag, Berlin, 2005), p. 185. 
2« H. D. Chen and Z. Nussinov, J. Phys. A 41, 075001 (2008); 

D. H. Lee, G.-M. Zhang, and T. Xiang, Phys. Rev. lett. 

99, 196805 (2007). 

2*^ S. Mondal, K. Sengupta and D. Sen, arXiv:0808.1175' 
^° A. Kopp and S. Chakravarty, Nature Physics 1, 53 (2005). 



